Sanorama Hie et al., 2019
GitHub
Tutorial external API
External external API tutorial
A fix to run scran pooling normalization computeSumFactors in current python environment.
import scanpy as sc
import scanpy.external as sce
import numpy as np
import pandas as pd
import os
# Working directory
os.chdir('/research/peer/fdeckert/FD20200109SPLENO')
# rpy2
os.environ['R_HOME'] = '/home/fdeckert/bin/miniconda3/envs/p.3.8.12-FD20200109SPLENO/lib/R'
# Plotting
import rpy2.robjects as robjects
color_load = robjects.r.source('plotting_global.R')
color = dict()
for i in range(len(color_load[0])):
color[color_load[0].names[i]] = {key : color_load[0][i].rx2(key)[0] for key in color_load[0][i].names}
sc.set_figure_params(figsize=(5, 5))
# Scanpy
n_comps=100 #PCA
dimred=n_comps #neighbors
n_neighbors=50 #neighbors
adata = sc.read_h5ad('data/object/so_sct.h5ad')
adata = adata.raw.to_adata()
def set_color(categories):
categories = [x for x in categories if x in list(adata.obs.columns)]
for category in categories:
adata.obs[category] = pd.Series(adata.obs[category], dtype='category')
keys = list(color[category].keys())
keys = [x for x in keys if x in list(adata.obs[category])]
adata.obs[category] = adata.obs[category].cat.reorder_categories(keys)
adata.uns[category+'_colors'] = np.array([color[category].get(key) for key in keys], dtype=object)
# Set colors
set_color(list(color.keys()))
sc.pp.normalize_total(adata)
sc.pp.log1p(adata)
sc.pp.scale(adata)
sc.tl.pca(adata, n_comps=n_comps)
adata.obs['integrate'] = adata.obs['treatment'].astype(str)+adata.obs['sample_rep'].astype(str)
sce.pp.scanorama_integrate(adata, key='integrate', basis='X_pca')
[[0. 0.74148021 0.45066851 0.27430691] [0. 0. 0.41067126 0.72495697] [0. 0. 0. 0.68332683] [0. 0. 0. 0. ]] Processing datasets NaClRep2 <=> NaClRep1 Processing datasets NaClRep1 <=> CpGRep1 Processing datasets CpGRep2 <=> CpGRep1 Processing datasets NaClRep2 <=> CpGRep2 Processing datasets NaClRep1 <=> CpGRep2 Processing datasets NaClRep2 <=> CpGRep1
# # Dimensional reduction and clustering
sc.pp.neighbors(adata, n_neighbors=n_neighbors, n_pcs=dimred, use_rep='X_scanorama')
sc.tl.leiden(adata, resolution=1)
sc.tl.louvain(adata, resolution=1)
sc.tl.umap(adata, min_dist=0.3)
# Plot
sc.pl.umap(adata, color=['louvain', 'leiden', 'tissue', 'treatment', 'label_fine_haemosphere', 'sample_rep', 'cc_phase_class', 'pHb_RNA', 'pRb_RNA', 'pMt_RNA', 'nCount_RNA', 'nFeature_RNA'], wspace=0.5, ncols=3)
adata = sc.read_h5ad('data/object/so_sct.h5ad')
hvg_8000 = list(adata.uns['hvg_int_8000'])
adata = adata.raw.to_adata()
def set_color(categories):
categories = [x for x in categories if x in list(adata.obs.columns)]
for category in categories:
adata.obs[category] = pd.Series(adata.obs[category], dtype='category')
keys = list(color[category].keys())
keys = [x for x in keys if x in list(adata.obs[category])]
adata.obs[category] = adata.obs[category].cat.reorder_categories(keys)
adata.uns[category+'_colors'] = np.array([color[category].get(key) for key in keys], dtype=object)
# Set colors
set_color(list(color.keys()))
sc.pp.normalize_total(adata)
sc.pp.log1p(adata)
sc.pp.scale(adata)
adata = adata[:,hvg_8000]
sc.tl.pca(adata, n_comps=n_comps)
adata.obs['integrate'] = adata.obs['treatment'].astype(str)+adata.obs['sample_rep'].astype(str)
sce.pp.scanorama_integrate(adata, key='integrate', basis='X_pca')
[[0. 0.77280551 0.43706777 0.30671613] [0. 0. 0.42271945 0.72908778] [0. 0. 0. 0.71339321] [0. 0. 0. 0. ]] Processing datasets NaClRep2 <=> NaClRep1 Processing datasets NaClRep1 <=> CpGRep1 Processing datasets CpGRep2 <=> CpGRep1 Processing datasets NaClRep2 <=> CpGRep2 Processing datasets NaClRep1 <=> CpGRep2 Processing datasets NaClRep2 <=> CpGRep1
# # Dimensional reduction and clustering
sc.pp.neighbors(adata, n_neighbors=n_neighbors, n_pcs=dimred, use_rep='X_scanorama')
sc.tl.leiden(adata, resolution=1)
sc.tl.louvain(adata, resolution=1)
sc.tl.umap(adata, min_dist=0.3)
# Plot
sc.pl.umap(adata, color=['louvain', 'leiden', 'tissue', 'treatment', 'label_fine_haemosphere', 'sample_rep', 'cc_phase_class', 'pHb_RNA', 'pRb_RNA', 'pMt_RNA', 'nCount_RNA', 'nFeature_RNA'], wspace=0.5, ncols=3)
adata = sc.read_h5ad('data/object/so_sct.h5ad')
hvg_6000 = list(adata.uns['hvg_int_6000'])
adata = adata.raw.to_adata()
def set_color(categories):
categories = [x for x in categories if x in list(adata.obs.columns)]
for category in categories:
adata.obs[category] = pd.Series(adata.obs[category], dtype='category')
keys = list(color[category].keys())
keys = [x for x in keys if x in list(adata.obs[category])]
adata.obs[category] = adata.obs[category].cat.reorder_categories(keys)
adata.uns[category+'_colors'] = np.array([color[category].get(key) for key in keys], dtype=object)
# Set colors
set_color(list(color.keys()))
sc.pp.normalize_total(adata)
sc.pp.log1p(adata)
sc.pp.scale(adata)
adata = adata[:,hvg_6000]
sc.tl.pca(adata, n_comps=n_comps)
adata.obs['integrate'] = adata.obs['treatment'].astype(str)+adata.obs['sample_rep'].astype(str)
sce.pp.scanorama_integrate(adata, key='integrate', basis='X_pca')
[[0. 0.76557659 0.46795758 0.33541585] [0. 0. 0.4454389 0.7555938 ] [0. 0. 0. 0.72159313] [0. 0. 0. 0. ]] Processing datasets NaClRep2 <=> NaClRep1 Processing datasets NaClRep1 <=> CpGRep1 Processing datasets CpGRep2 <=> CpGRep1 Processing datasets NaClRep2 <=> CpGRep2 Processing datasets NaClRep1 <=> CpGRep2 Processing datasets NaClRep2 <=> CpGRep1
# # Dimensional reduction and clustering
sc.pp.neighbors(adata, n_neighbors=n_neighbors, n_pcs=dimred, use_rep='X_scanorama')
sc.tl.leiden(adata, resolution=1)
sc.tl.louvain(adata, resolution=1)
sc.tl.umap(adata, min_dist=0.3)
# Plot
sc.pl.umap(adata, color=['louvain', 'leiden', 'tissue', 'treatment', 'label_fine_haemosphere', 'sample_rep', 'cc_phase_class', 'pHb_RNA', 'pRb_RNA', 'pMt_RNA', 'nCount_RNA', 'nFeature_RNA'], wspace=0.5, ncols=3)
adata = sc.read_h5ad('data/object/so_sct.h5ad')
hvg_4000 = list(adata.uns['hvg_int_4000'])
adata = adata.raw.to_adata()
def set_color(categories):
categories = [x for x in categories if x in list(adata.obs.columns)]
for category in categories:
adata.obs[category] = pd.Series(adata.obs[category], dtype='category')
keys = list(color[category].keys())
keys = [x for x in keys if x in list(adata.obs[category])]
adata.obs[category] = adata.obs[category].cat.reorder_categories(keys)
adata.uns[category+'_colors'] = np.array([color[category].get(key) for key in keys], dtype=object)
# Set colors
set_color(list(color.keys()))
sc.pp.normalize_total(adata)
sc.pp.log1p(adata)
sc.pp.scale(adata)
adata = adata[:,hvg_4000]
sc.tl.pca(adata, n_comps=n_comps)
adata.obs['integrate'] = adata.obs['treatment'].astype(str)+adata.obs['sample_rep'].astype(str)
sce.pp.scanorama_integrate(adata, key='integrate', basis='X_pca')
[[0. 0.76592083 0.44744122 0.35259664] [0. 0. 0.46678141 0.77487091] [0. 0. 0. 0.71261226] [0. 0. 0. 0. ]] Processing datasets NaClRep1 <=> CpGRep1 Processing datasets NaClRep2 <=> NaClRep1 Processing datasets CpGRep2 <=> CpGRep1 Processing datasets NaClRep1 <=> CpGRep2 Processing datasets NaClRep2 <=> CpGRep2 Processing datasets NaClRep2 <=> CpGRep1
# # Dimensional reduction and clustering
sc.pp.neighbors(adata, n_neighbors=n_neighbors, n_pcs=dimred, use_rep='X_scanorama')
sc.tl.leiden(adata, resolution=1)
sc.tl.louvain(adata, resolution=1)
sc.tl.umap(adata, min_dist=0.3)
# Plot
sc.pl.umap(adata, color=['louvain', 'leiden', 'tissue', 'treatment', 'label_fine_haemosphere', 'sample_rep', 'cc_phase_class', 'pHb_RNA', 'pRb_RNA', 'pMt_RNA', 'nCount_RNA', 'nFeature_RNA'], wspace=0.5, ncols=3)
adata = sc.read_h5ad('data/object/so_sct.h5ad')
hvg_2000 = list(adata.uns['hvg_int_2000'])
adata = adata.raw.to_adata()
def set_color(categories):
categories = [x for x in categories if x in list(adata.obs.columns)]
for category in categories:
adata.obs[category] = pd.Series(adata.obs[category], dtype='category')
keys = list(color[category].keys())
keys = [x for x in keys if x in list(adata.obs[category])]
adata.obs[category] = adata.obs[category].cat.reorder_categories(keys)
adata.uns[category+'_colors'] = np.array([color[category].get(key) for key in keys], dtype=object)
# Set colors
set_color(list(color.keys()))
sc.pp.normalize_total(adata)
sc.pp.log1p(adata)
sc.pp.scale(adata)
adata = adata[:,hvg_2000]
sc.tl.pca(adata, n_comps=n_comps)
adata.obs['integrate'] = adata.obs['treatment'].astype(str)+adata.obs['sample_rep'].astype(str)
sce.pp.scanorama_integrate(adata, key='integrate', basis='X_pca')
[[0. 0.75215146 0.44190871 0.34923928] [0. 0. 0.51669535 0.76351119] [0. 0. 0. 0.76513081] [0. 0. 0. 0. ]] Processing datasets CpGRep2 <=> CpGRep1 Processing datasets NaClRep1 <=> CpGRep1 Processing datasets NaClRep2 <=> NaClRep1 Processing datasets NaClRep1 <=> CpGRep2 Processing datasets NaClRep2 <=> CpGRep2 Processing datasets NaClRep2 <=> CpGRep1
# # Dimensional reduction and clustering
sc.pp.neighbors(adata, n_neighbors=n_neighbors, n_pcs=dimred, use_rep='X_scanorama')
sc.tl.leiden(adata, resolution=1)
sc.tl.louvain(adata, resolution=1)
sc.tl.umap(adata, min_dist=0.3)
# Plot
sc.pl.umap(adata, color=['louvain', 'leiden', 'tissue', 'treatment', 'label_fine_haemosphere', 'sample_rep', 'cc_phase_class', 'pHb_RNA', 'pRb_RNA', 'pMt_RNA', 'nCount_RNA', 'nFeature_RNA'], wspace=0.5, ncols=3)